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ABSTRACT 

X-ray sources with very few counts can be identified with low-noise X-ray 
detectors such as the Advanced CCD Imaging Spectrometer onboard the Chandra 
X-ray Observatory. These sources are often too faint for parametric spectral 
modeling using well-established methods such as fitting with XSPEC. We discuss 
the estimation of apparent and intrinsic broad-band X-ray fluxes and soft X-ray 
absorption from gas along the line-of-sight to these sources, using nonparametric 
methods. Apparent flux is estimated from the ratio of the source count rate 
to the instrumental effective area averaged over the chosen band. Absorption, 
intrinsic flux, and errors on these quantities are estimated from comparison of 
source photometric quantities with those of high signal-to-noise spectra that were 
simulated using spectral models characteristic of the class of astrophysical sources 
under study. 

The concept of this method is similar to the long-standing use of color- 
magnitude diagrams in optical and infrared astronomy, with X-ray median energy 
replacing color index and X-ray source counts replacing magnitude. Our nonpara- 
metric method is tested against the apparent spectra of ~ 2000 faint sources in 
the Chandra observation of the rich young stellar cluster in the M 17 HII region. 
We show that the intrinsic X-ray properties can be determined with little bias 
and reasonable accuracy using these observable photometric quantities without 
employing often uncertain and time-consuming methods of non-linear parametric 
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spectral modeling. Our method is calibrated for thermal spectra characteristic 
of stars in young stellar clusters, but recalibration should be possible for some 
other classes of faint X-ray sources such as extragalactic active galactic nuclei. 

Subject headings: methods: data analysis - methods: statistical - open clusters 
and associations: individual (M17) - techniques: photometric - X-rays: general - 
X-rays: stars 



1. INTRODUCTION 



Many thousands of X-ray sources are discovered using X-ray cameras with modest spec- 
tral resolution, such as the CCD detectors on the ASCA, Chandra, and X MM- Newton 
space observatories. After detecting and locating the source, the observer often seeks to es- 
timate intrinsic properties of the source, particularly its flux in a broad band of interest and 
spectral characteristics such as power law index, thermal plasma temperature, or column 
density of line-of-sight gas absorbing soft X-rays. Low-resolution source and background 
spectra are extracted from the CCD data, and x 2 ot maximum likelihood statistical nonlin- 
ear regression of these data with multi-parameter spectral models is performed with codes 
such as XSPEC0, Sherpjfl, or MIDAsjE 

We have been particularly concerned with the study of large populations of faint X- 
ray sources as sociated with pre-ma i n sequence (PMS) stars residing in or near molecu 



lar cl ouds (e.g. iFeigelson et al.l 120021 : iGetman et al.ll2005l : iTownsley et al.ll2006l : iBroos et al. 



20071 ). These and other studies of X-ray-bright young stars use a spectral model family 
involving one- or two-tem perature thermal p lasmas with a particular pattern of non-solar 



elemental abundances (e.g. lMaggio et al.ll2007l ). Each star is subject to a different absorption 



due to its chance location in the molecular cloud and, for the youngest systems, its local pro- 
tostellar envelope. However, when faint sources are considered, there can be too few counts 
to obtain reliable spectral fits using regression methods; a wide range of models can be fit to 
the same distribution of source photon energies. While apparent broad-band fluxes can still 
be estimated, the role of absorption, the intrinsic fluxes, and true plasma temperature(s) are 
often poorly constrained from spectral modeling of faint sources. Also, parametric regression 
analyses are often inhomogeneous in practice; two-temperature plasma models may be fit 



1 http : //heasarc . gsf c . nasa.gov/docs/xanadu/xspec/ 
2 http : // cxc . harvard . edu/ sherpa/ 

http : //www. eso . org/sci/data-processing/sof tware/esomidas/ 
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to lightly-absorbed, bright (say, >200 count) sources while fainter and/or highly-absorbed 
sources may be fit with one-temperature plasma models. There is thus a motivation to con- 
sider nonparametric approaches to interpreting X-ray CCD energy distributions which can 
be applied in a uniform fashion to bright and faint sources. 

In this work, we examine the nonparametric inference of intrinsic X-ray luminosities 
and absorbing column densities which can be efficiently obtained for large populations of 
both bright and faint X-ray CCD sources. The concept of the method is similar to the long- 
standing use of color-magnitude diagrams in optical and infrared astronomy. For PMS stars, 
a photometric magnitude and at least one suitable color index, such as the J vs. J — H color 
magnitude diagram, are used in conjunction with distance, age estimates, and stellar interior 
evolutionary models to estimate stellar masses, bolometric luminosities, and line-of-sight ab- 
sorption. In our proposed method, photometric X-ray quantities such as apparent hard-band 
background-corrected flux (analogous to the J magnitude) and background-corrected total- 
band median energy (analogous to the J — H color) are used to estimate absorption-corrected 
X-ray fluxes, and line-of-sight absorbing column densities. We apply our nonp arametric es 



timates to a large sample of faint X-ray sources in the M 17 star-forming region (iBroos et al. 



20071 . and Townsley et al., in preparation). Intrinsic spectra for this class of stars are cali- 
brated to the high signal-to -noise Chandra ACI S-I spectra from the sensitive Chandra Orion 
Ultradeep Project (COUP; iGetman et al.ll2005l ). Reliability of our nonparametric methods, 
including both systematic bias and statistical uncertainty, are evaluated using simulated 
X-ray sources with known properties. 

The technical concept of avoiding nonlinear par ametric modeling, e specially in cases of 
faint sources, is not new. In Chandra ACIS studies, iGagne et al.l (120041 ) used the extracted 
X-ray mean energy and its comparison with the simulated spectra of one-temperature models 
to infer some stellar properties. Our method differs in the use of more physically realistic 
spectral models a nd the use of the m edian instead of mean energy to be statistically robust 
against outliers. iHong et al.l (120041 ) also comment that median energies can give effective 
spectral estimators in X-ray CCD spectroscopy, but they use the median in conjunction 
with quartiles which can be statistically less reliable for very faint sources. For example, 
the 25%-50%-75% energy quartiles are readily estimated for a source with 9 photons (the 
energies of the 3 rd , 5 th , and 7 th photons), but are not readily estimated if 8 or 10 photons 
are present. 

This study is closely linked to the ACIS Extract softw are package for Chandra ACIS 
data analysis which is described by iBroos et al.l (120021 . 120091 ) . We use ACIS Extract version 



http : //www. astro .psu . edu/xray/docs/TARA/ae_users_guide . html 
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2009-01-27 in the analysis here. 



Methods 



2.1. Spectral Models for PMS Stars 

An initial step requires the definition of the intrinsic spectral shape of the population 
of X-ray sources under study; this may be flux- dependent. This shape may be obtained 
from previous spectral analysis of brighter sources of the same class, or from astrophysical 
models. In our application to young stars, we use flux-dependent thermal plasma models. 
In other applications, such as active galactic nuclei in extragalactic surveys or X-ray binaries 
in nearby galaxy studies, different spectral shapes would be used (j2J. 

We use the COUP spectral results as calibrators of ACIS observations of star populations 
in more distant star formation regions. Spectral templates derived from bright COUP sources 
will be used to translate apparent photometric properties of weak PMS stars into their 
intrinsic properties. In §2.21 and §2.31 we will concentrate on two observable properties 
of each source: the apparent hard-band X-ray flux, F^, and the X-ray total-band median 
energy of background-corrected extracted counts, MedE t . Here the hard-band covers the 
2.0 — 8.0 keV energy range and the total-band covers the 0.5 — 8.0 keV energy range. MedE t 
and a nonparametric estimator of are auto matically calculated from the extracted source 
photon events by the A CIS Extract package (IBroos et al.l 120021 . 120091 ) . Later in the study, 
we consider total-band fluxes and absorption-corrected fluxes. 

We use the COUP spectral fits to establish a p arametric spectral mo del family as a 
basis for estimating broad-band fluxes of PMS stars. iPreibisch et al.l (120051 ) find that two- 
temperature thermal plasma models provide good fits to high signal-to-noise ACIS data of 
~ 500 COUP PMS stars. The temperature of the hot plasma component and emission 
measure ratio between hot and cool components increase with stellar surface flux, with 
the temperature of the cool plasma component remaining approximately constant at kT\ ~ 
0.8 keV (T ~ 9 MK). While it is recognized that the coronae of active stars intrinsically have 



a distribution of plasma temperatures (e.g. iGiidel et al.l 120071 ) . the two-temperature model 



is generally adequate for most ACIS-quality spectra. COUP PMS stars follow the general 
correlation between the temperatures of the hot an d cool plasma compon ents seen in main- 
sequence stars, but with much higher temperatures (IPreibisch et al.ll2005l ). The 0.8 keV cool 
plasma component likely describes a persistent coronal structure of PMS stars in the form 
of compact loops with high plasma density, while the ho t plasma component reflects the 
contribution of PMS flaring activity (IPreibisch et al.ll2005l . and references therein). 
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Figure [I] presents the dependencies of the hot plasma component, kT 2 , and the ratio 
of hot to cool component emission measures, EM 2 / EMi, on the absorption-corrected hard- 
band luminosit y , Lh. r , for COUP PMS stars. A similar diagram appears in Figure 11 of 



n,ci 

Preibisch et all (120051) . A few d ozen COUP super-hot stars with their kT 2 above 5 keV 



are omitted (IGetman et al.l 120081 ). The square symbols in the figure and the parameters in 
Table [1] define a set of eleven model families we adopt for the class of PMS stars. These are 
median values with linear extrapolations at the ends of the range]. The dashed lines 
around the median values encompass ~ 75% of COUP stars and indicate the uncertainty 
in choosing an X-ray model among the model families we consider. This information will 
be used below to estimate the sensitivity of derived stellar properties on assumed spectral 
models. 

We note that the XEST project reports little correlation of plasma temperature with 
X-ray luminosity in Class II stars (PMS stars with accretion disks), and higher plasma 
temperatures in Clas s II compared to Class III stars (PMS stars with weak or absent disks) 
( jTelleschi et al.ll2007l ). The COUP sample similarly shows a larger fraction of active accretors 
lying outside of the band of the X-ray model uncertainty (FigureEJ). Thus the spectral models 
used here best describe Class III stars, which generally dominate X-ray-selected samples of 
young rich X-ray stellar populations. The more general point is that care must be taken 
to match the spectral model to the underlying population when photometric quantities are 
examined. 

We simula te 11 COUP model families from Tabled] using the fakeit command in XSPEC 
( jArnaudl 119961 ). Each model family is characterized by values of kT 2 and EM 2 f EM\ with 
fcTi = 0.8 keV and is simulated using the WABS x (MEKAL+MEKAL) XSPEC model 
family with 0.3 solar elemental abundances characteristic of young starfl For each model 
family, we generate a large number of high signal-to-noise Chandra ACIS-I spectra over a 
grid of absorption column densities within the range of 20 < log(iV#) < 24 cm -2 . Simulated 
spectra are then passed through the ACIS Extract code for photometric analysis including 



5 The spectral trend is compatible with th e L x ex T 4 ;. 3 ~ 6 7 re lation (jTelleschi et al.ll2007h seen in the XMM- 
Newton Extended Survey of Taurus (XEST; lGudel et al.ll2007l ). where T av is the mean electron temperature 



weighted by the differential emission measure distribution and Lx is the intrinsic X-ray luminosity derived 
in the 0.3 — 10 keV energy band. 



6 This simulation was chosen to be compatible with the COUP study (jGetman et al.l I2005T ) where 
bright lightly ab sorbed sour ces were modeled with two-temperature optically thin thermal plasma 



MEKAL models (Mewe 1991), assuming a uniform density plasm a with 0.3 times solar elemental abun- 
dances jlm anishi. Kovama. fc Tsuboil 120011 iFeigelson et al.l 12002). Solar abundances were taken from 
Anders fc Grevesse Jl989l). X-ray a"bsorption was modeled using the WABS model of atomic cross sections 
of iMorrison fc McCammonl |l983h . 
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calculation of source net counts, median energies, and apparent fluxes in several energy 
bands. 



2.2. Apparent Photometric Fluxes from ACIS Extract 

The ACIS Extract (AE) package has been used in do zens of Chandra- A CIS studies, 
particularly projects involving analysis of many faint sources faroos et al.ll200gh . It computes 



several photometric quantities of interest here from the extracted source events. NCt and 
NCh (in counts) are the total-band and hard-band net counts after subtraction of local 
background. MedE t and MedE^ (in keV) are the total-band and hard-band apparent X-ray 
median energies after background subtraction. The program estimates photometric photon 
fluxes (in photons cm~ 2 s _1 ) in the total and hard bands using 

Ft, pho t = NC t / <ARF t > /Exp (1) 
F h)Phat = NC h / <ARF h > /Exp. (2) 

Here < ARE > is the mean value of a source's Auxiliary Response File (ARF), the effective 
area averaged over the appropriate band, and Exp is the exposure time at the source location 
on the detectoi0. We estimate the apparent energy fluxes F t and Fh (in ergs cm -2 s _1 ) as 

F t = 1.602 x 1(T 9 x F t , phot x MedE t (3) 
F h = 1.602 x 10- 9 x F Kphot x MedE h . (4) 

where MedE are in keV and the constant comes from the keV to ergs conversion. Luminosi- 
ties L t and Lh are derived from these fluxes by multiplying them by AirD 2 where D is the 
object's distance, in cm. 

In various analyses below, we will consider fluxes that have been corrected for soft X- 
ray absorption scaled to the line-of-sight column density Nh which we derive from MedE t . 
These intrinsic absorption-corrected flux estimates are labeled F tc and Fh c (in ergs cm -2 s _1 ) 
for the total and hard bands, respectively. 

In oder to assess the fidelity of these flux estimates, Figure [3] compares ACIS Extract 
photometric fluxes from equations (3)- (4) for simulated data with apparent fluxes derived 



7 In the ACIS Extract documentation, F t , p hot is designated F2, while F\ represents another estima- 
tor of incident flux where effective areas are measured in narrow ene rgy bands. The Fi-type fluxes are 
preferred for bright sources but are less stable at low count rates (jBroos et al.l 120021 120091 ) . Fi-type 
fluxes are used in the Chandra Source Catalog where they are called 'Aperture Source Energy Fluxes' 
(http : //asc . harvard . edu/ csc/columns/f luxes . html) . 
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from the flux integration of the input models (discussed in §2.11) evaluated on the energy bins 
set by the Chandra response matrix. Discrepancies of 10 — 20% in flux are seen for typical 
spectra, and can range up to ~ 40% for extremely soft or hard spectra. These biases can be 
explained by inaccuracies in the computation of fluxes using equations (l)-(2) that are based 
on effective areas averaged over the broad band assuming a flat incident spectrum. This 
approximation is not accurate for realistic spectral shapes. ACIS Extract iq-type incident 
fluxes, where effective areas are measured in narrow energy bands, do not exhibit these 
biases, but they suffer from large Poisson errors for weak sources. For example, a single 
photon present at energies E > 6 keV, where the Chandra mirror response is poor, can lead 
to a large spurious jump in inferred flux. Since our work is oriented towards study of faint 
sources, we use the more stable F 2 -type photometric fluxes of equations (l)-(2). The curves 
in Figure [3] can then be used to reduce the systematic biases in flux extimates for extremely 
soft and hard spectra (see also Table [2]). 



2.3. Association Between Apparent Median Energy and Inferred Absorption 

F t and Fh are not reliable estimators of the intrinsic source flux if absorption is present. 
For bright sources, a typical procedure is to estimate line-of-sight column density Nh via a 
parametric fit of the apparent photon energy distribution using XSPEC or similar code, and 
calculate an intrinsic unabsorbed flux based on the other spectral model parameters. We 
now seek an analogous procedure based on the apparent source median energy, MedE t . 

Figure|U(see also Table[3]) shows the simulated X-ray column density (N H ) as a function 
of the X-ray source's median energy in t he total band {M e dEt) . The results of our simulations 



confirm the previous COUP finding of iFeigelson et al.l (120051 ) that the median energy can 



be effectively used as a surrogate for absorption column density for MedE t > 1.5 — 1.7 keV 
corresponding to log Nh > 22.0 cm -2 . The 0.5 keV low energy limit of the ACIS instrument 
leads to a degeneracy in absorption values for MedE t ~ 1.0 — 1.3 keV. The MedE t — log Nh 
relationship can be used reliably to infer either absorbing column density or source plasma 
temperature if Nh is known in advance. If both are unknown, the inference has greater error 

(mm- 

Figure [5] (see also Table H|) shows the ratio of intrinsic (absorption-corrected) flux to 
apparent flux against apparent median energy in the total band for the simulated spectra. 
These are flux integrations of the input models (not re-fitted simulated spectra) evaluated 
on the energy bins set by the Chandra response matrix. We see that flux ratios for the total 
band fluxes, F tc /F t , are often factors of ~ 10, and can reach factors of ~ 30, for heavily 
absorbed sources. In the hard band, the flux ratios F^ c j F^ are much smaller, often less than 
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a factor of 2 and always less than a factor of 10. Fh c is thus a much more stable estimate of 
source flux than F tc when absorption is present. 



2.4. Implementation 

As described above, the inferred relation between the intrinsic spectral shape and the X- 
ray luminosity ( §2.11 and Figures [1] and [2]) allowed us to explicitly define and simulate X-ray 
models of PMS stars. Analysis of the simulated data calibrates relationships between spectral 
models and observable quantities ( §2.2l and l2.3[ and Figures El HI EI). We have shown that the 
source counts and median energy, which are readily measured for faint sources, often serve as 
surrogates for source flux and absorption. Similar observable quantities have been used for 
many years in color-magnitude diagrams for characterizing sources in optical-infrared (OIR) 
astronomy. 

Other classes of sources, different from young stars, may have their own relationships 
between intrinsic and/or observable quantities that allow simulations of explicit X-ray models 
O- Once model simulations are performed, calibration data mapping simulated observable 
to simulated intrinsic properties can be built. In the case of PMS stars we produce four 
calibration products: Tables [2j [3j and H] correspond to Figures IHE discussed above, and 
Table is used for error analysis ( §2.4.21) . 

Once simulated calibration data are built, the measured photometric properties of ob- 
served X-ray sources are compared with the calibration data to obtain intrinsic properties 
and errors on these properties. The procedure is shown in Figures El and [7] and is pre- 
sented in detail in this section. The IDL implementation of the procedure can be found 



at http://www.astro.psu.edu/users/gkosta/XPHOT/. Researchers working with other 



classes of sources can use the IDL code after constructing the four calibration tables similar 
to Tables [2]l5] appropriate for the spectral shapes of their source class. 



2.4-1- Translating Apparent to Intrinsic Properties 

In Step 1 of the procedure shown in Figure El source photometry properties are obtained 
by running ACIS Extract. Here we provide details on Steps 2-4 which estimate the source's 
absorption and intrinsic broad-band fluxes. 

In Step 2, one of the 11 simulated PMS spectral model families listed in Table [T] is 
assigned to each individual source of interest based on that source's hard-band apparent 
luminosity and full-band median energy. For lightly obscured stars (MedE t < 2 keV), the 
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ratio of intrinsic to apparent flux in the hard-band is < 0.2 dex for all spectral models 
(Figure [5b), so that apparent hard-band luminosity L h only slightly underestimates the 
intrinsic hard-band source luminosity Lh c - The < 0.2 dex luminosity difference is much less 
than the quantization of 0.5 dex on Lh c used to create the COUP template models (Figure 
[Hand Tabled]). Thus, for a lightly obscured PMS star, the apparent luminosity Lh can be 
used directly to assign a specific spectral model family ("best model family") from those 
listed in Table [H 

For a heavily obscured PMS star (MedE t > 2 keV), the source's median energy MedE t 
is compared to the calibration of intrinsic to apparent flux ratios (Figure [5b and Table H]) 
to choose 11 candidate F hc / F h values corresponding to 11 template model families. Then, 
the source's apparent hard-band luminosity L h is used to obtain 11 candidate intrinsic hard- 
band luminosities, Lh c ,candidate- The "best model family" is the one with the closest match 
between the Lh c ,candidate and the "nominal" intrinsic hard-band luminosity listed in Table [TJ 
This procedure is analogous to the reddening correction in OIR color-magnitude diagrams. 

In Step 3, photometric fluxes are de-biased using calibration curves (Figure [3] and Table 
[2]) for the source's spectral model. The absorbing column density is then estimated using 
the source's median energy (Figure 0] and Table [3]). Step 4 uses the resulting Nh estimate 
to obtain absorption-corrected fluxes (Figure [5] and Tabled]). 

To visually represent the transformation of apparent to intrinsic stellar properties in 
our method, Figure [S] shows the X-ray "color-magnitude diagram" derived from the simu- 
lated PMS star data. In this diagram, the apparent hard-band luminosity (Lh), apparent 
total-band median energy (MedEt), absorbing column density (Nh), and intrinsic hard-band 
luminosity (Lhc) are analogues to OIR magnitude, color, extinction, and bolometric lumi- 
nosity, respectively. We note again the degeneracy in low values of absorbing column density 
at MedE t ~ 1 — 1.3 keV due to the 0.5 keV limit of the Chandra- AC1S instrument. 



2.4-2. Error Analysis 

We estimate uncertainties in source properties derived using our photometric procedures 
including both small- N statistical errors and systematic uncertainties in our knowledge of the 
X-ray spectral model. Systematic errors associated with the astronomical source properties, 
such as distance or X-ray variability, are not considered here. We also omit systematic errors 
of instrumental origin. Errors in the ARF, which is based on calibration of the Chandra 
mirrors and ACIS detector, are anticipated to be sma ll and become com parable to statistical 
uncertainties only for sources with > 10, 000 counts (jDrake et al. 2006 ). The procedures for 
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the error analysis of photometric quantities described here are outlined in Figure [7J 

First, we estimate statistical errors on the X-ray median energy statistic using simulated 
data. Passing simulated model spectra ( §2.ip through the MARX mirror-detector simula- 
we slice the resulting X-ray event lists into thousands of subsamples of different sizes. 
Distributions of the median energy statistic are then obtained and corresponding standard 
errors, AMedE, encompassing 68% of the values are recorded. The resulting standard errors 
for the total-band and hard-band median energies for the 11 PMS spectral model families 
are shown in Figure M and tabulated in Table 

Figure [9] shows that median energy errors range from ~ 5% for sources with 100 counts 
to ~ 15% for sources with 10 counts. These values are larger than expected from a standard 
Gaussian distribution, although the errors do scale roughly with \/N. At very low (MedE t < 
2 keV) and very high (MedE t > 5.5 keV) median energy values, errors are reduced due to 
the energy limits of the Chandra instrumental response. The spread in AMedE t for a fixed 
source count and MedE t value shows that sources with harder intrinsic spectra have larger 
median energy uncertainties than intrinsically softer sources^]. 

Again we recall that the results here are specific to the Chandra ACIS mirror-detector 
combination and PMS stellar spectral models. Researchers working with data from other 
X-ray telescopes or source classes are advised to produce new error estimates for the median 
energy statistic. However, we expect that the qualitative findings will be similar for other 
systems. 

To estimate the effects of these median energy uncertainties (Figure [9] or Table EJ), we 
propagate them along with uncertainty on s ource net counts to errors on apparent fluxes 



(equations 1-4) using standard methods (e.g. iBevington fc Robinson! I1992T ). This is Step 5 
of the procedure shown in Figure [7J 

Step 6 of the procedure propagates the median energy errors through the relations 



s http : //space . mit . edu/CXC/MARX/ 

9 We also checked our error results using a different method which directly employs generated spectra 
without the need of MARX simulations. Here we generate a large number of simulated spectra with the 
desired number of counts using random deviates from the known energy distributions for each spectral model 
( Press et al 1 ll992h . The resulting median energy distributions and standard error values are indistinguishable 



from those obtained from the MARX simulations. 

10 When a low-temperature and a high-temperature source have the same observed median energy, the 
low-temperature source has to be more highly absorbed (Figure 0}. That low-energy Nji cut-off and the 
intrinsic lack of power at high energies give the low-temperature source a narrower observed spectrum, and 
a correspondingly more accurate estimator for apparent median energy. 
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between median energy and absorbing column density (Figure H] and Table [3]) to estimate 
both statistical and systematic errors in log Nh- This is illustrated in the insert diagram 
of Figure HI Statistical errors on log N H are derived by propagating errors on MedE t using 
the curve for the "best model family" ( §2.4.1)) . Systematic errors on logiV# are estimated 
by propagating the value of the MedE t itself to other possible spectral models; we use the 
"second nearest neighbor" as a plausible range for model uncertainty (see the dashed lines 
in Figure [1] which are discussed in §2.ip . 

Finally, we estimate statistical and systematic errors on the absorption-corrected hard- 
band (Fhc) and full-band (F tc ) fluxes by propagating errors already obtained for MedE t , Fh, 
and F t . This propagation uses the absorption correction calibrations in Figure [5] and Table 
HI As with log Nh, statistical errors on intrinsic source fluxes are estimated using the best 
spectral model, while systematic errors are estimated using the "second nearest neighbor" 
models. 



3. Application to M17 X-ray Sources 



We now illustrate our analysis methods using the Chandra ACIS-I X-ray data of PMS 
candidate stars in the M17 region. M17, one of the brightest HII regions in the sky, is 
ionized by OB stars in the rich young stellar cluster NGC 6618 which lies on the edge of a 
massive molecular cloud. The X-ray stellar population is composed of light to moderately- 
absorbed PMS stars of the central NGC 6618 cluster supplemented by heavily-absorbed stars 
embedded in or obscured by the surrounding cloud. A study of the X-ray popul ation from 
a 40 ks Chandra ACIS exposure of the region is described by iBroos et al.l (120071 ); 886 faint 
X-ray sources were found in this short exposure. 

A new ~ 300 ks Chandra exposure of this field has been obtained, and the combined 
dataset yields ~ 2000 X-ray sources (Townsley et al., in preparation). This dataset is chosen 
for its large sample and wide range of source properties: 5 < NC t < 1000 source counts and 
1 < MedE t < 5 keV median ene rgies. The Chandra data have been analyzed following ACIS 
Extract procedures described in lBroos et al.l (120021 . 120091 ). Spectral fits have been performed 
with XSPEC on sources with net counts down to ~ 10 counts; best-fit models were found by 
the maximum likelihood method"] assuming a one-temperature TBABS x VAPEC model 



11 The C-statistic was used to fit ungrouped spectra for more than 2000 M17 sources. The back- 
ground spectra available for most sources contained ~ 100 counts. Within XSPEC these moderate- 
quality backgrounds were modeled with a simple continuous piecewise-linear function with 10 ap- 
proximately evenly-spaced vertices (see the discussion of the cplinear model in the AE manual, 
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( I Wilms et al.l I200CH : ISmith et al.l 1200 ll ) with the elemental abundances fixed at the value s 
typical for PMS or extremely active zero-age main-sequence stars from iGiidel et al.l (120071 ). 
These are slightly different spectral models to those used in the COUP study, but the small 
change in spectral shape is likely to have negligible effect on the inferred absorption and flux 
estimates. 

The Ml 7 data thus give a large sample of PMS X-ray sources for which spectral and 
photometric properties are available. The absorption and flux estimates derived with the 
parametric and nonparametric methods can thus be compared. 



3.1. X-ray column densities towards M17 sources 

Figure [TOl compares the nonparametric and parametric estimation of absorption for M17 
sources. Panel (a) shows the MedE t obtained nonparametrically plotted against the column 
density log Nh obtained from XSPEC. The data accurately follow the relationship shown in 
Figure H] with increased scatter for fainter sources, as expected. Panel (b) shows the ratio 
of Nh values estimated from MedE t as described in § 12.31 with those obtained from XSPEC 
fits, as a function of MedE t . For the MedE t > 1.7 keV range, the log Nh value inferred 
nonparametrically is systematically lower than that of the XSPEC fitting by 0.05 dex. This is 
a small effect, and we do not know if it represents a bias in our MedE t — log Nh conversion, a 
bias in XSPEC fitting, or is due to the difference in spectral models used in the simulations 
and M17 spectral fitting. For sources with MedE t < 1.7 keV, discrepancies between the 
nonparametric and parametric absorption estimates show much larger scatter and biases, 
up to a factor of ten. Neither method can accurately measure absorbing column density for 
faint soft sources, because Chandra ACIS-I spectra are often insensitive to differences in low 
column density and there is an ambiguity in spectral model fits of faint soft sources. 

We derive statistical and systematic errors on absorption column density following the 
procedure presented in §2.4.21 Figure [TT] shows that log Nh values inferred from our method 
have increasing uncertainties for soft sources with MedE t < 1.7 keV. Again, this arises from 
the weak dependency of apparent median energy on column density for softer sources with 



http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.pdf). In order to explore the range of al- 
ternative model fits the TBABS x VAPEC model with thawed parameters was fit to each source five times, 
using a set of initial parameter values that explore an appropriate region of parameter space. Two additional 
TBABS x VAPEC models representing "standard" stellar spectra (kT frozen at 0.86 and at 2.6 keV) were 
also fit. The best model derived from these seven fits was chosen via a combination of computer algorithm 
and human review. The full procedure requires days to be accomplished. 
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light obscuration (Figure H]). We recommend that tabulations of nonparametric properties 
of PMS stars with MedE t < 1.7 keV should report an upper limit of log(iV#) < 22.0 cm -2 
rather than specifying a very uncertain absorbing value. 

For stars with MedE t > 1.7 keV, Figure [TTa shows that our method estimates logiVfl- 
with a statistical accuracy better than ±0.3 dex, ±0.2 dex, ±0.1 dex for the NCt = 7 — 20, 
20 — 50, > 50 count strata, respectively. Panel (b) shows that systematic uncertainties due 
to inadequate knowledge of the precise X-ray spectral model applicable to each source add 
an additional uncertainty less than ±0.15 dex for all stars, and less than ±0.05 dex for 
sources with NCt ^ 100. These errors are comparable to the scatter between photometric 
and XSPEC estimates of logiV# above 22.0 cm -2 shown in Figure [TOb. In summary, the 
nonparametric procedure provides reliable estimates of Nh absorption within a factor of 
~ 1.5 for sources with > 50 counts, and within a factor of ~ 2 for sources with ~ 10 counts 
including both statistical and systematic spectral model uncertainties. 



3.2. Apparent X-ray fluxes of M17 sources 

Figure [T2h shows the ratio of the total-band apparent flux obtained using our nonpara- 
metric methods in § 12.21 and the flux obtained from XSPEC fitting for Ml 7 sources. The 
biases shown in Figure [3] have been removed from the nonparametric estimates. The mean 
value is 1.04 ± 0.10 indicating that only a small (< 4%) systematic flux difference between 
the nonparametric and parametric methods remains, and individual source differences are 
mostly below ±10%. In Figure [12b. similar distributions for hard-band apparent flux ratios 
are shown, again exhibiting small (~ 2%) bias. These small systematic differences between 
nonparametric and XSPEC results may become important in the > 100 — 300 counts regime, 
and parametric methods are preferred for strong sources. 

Following procedures outlined in § I2.4.2[ we derive statistical errors on apparent fluxes 
and plot them in Figure O Panel a shows that our method calculates apparent total-band 
source fluxes with statistical accuracy better than 60%, 50%, 30%, 20% for total-band net 
count strata of 7 — 10, 10 — 20, 20 — 50, > 50, respectively. Panel h shows similar errors for 
apparent hard-band source fluxes. These are slightly higher than Gaussian yN errors, and 
reach \fN levels when > 100 counts are present. 
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3.3. Absorption Corrected X-ray Fluxes of M17 Sources 

Figures [T4"lfT5l are similar to Figures [T2lfT3l but for broad-band intrinsic fluxes of M17 
sources corrected for absorption. Recall that very large corrections have been applied to high- 
MedE t sources (Figure[5]). Here we see considerably greater bias and scatter compared to the 
apparent fluxes; note the expanded vertical scale in Figure [T4"b. The ratio of nonparametric 
to XSPEC intrinsic total-band fluxes is 1.45 ± 0.74 (1.34 ± 0.19) for the 20 - 50 (> 300) 
total-band count strata, and the ratio of hard-band fluxes is 1.10 ± 0.35 (1.07 ± 0.04) for 
the 20 — 50 (> 300) hard-band count strata. Thus, even for high-count rates, the intrinsic 
photometric total-band (hard-band) fluxes are systematically ~ 35% (~ 7%) higher than 
XSPEC fluxefl 

For any broad band flux that includes the soft 0.5 — 2 keV band, the correction from ap- 
parent to intrinsic fluxes can be both very large (Figured) and very uncertain (Figure [15b), 
even for moderate absorption. Weaker sources are particularly vulnerable. We discourage 
use of the F tc and L tc quantities derived from photometry for scientific analysis. Although 
we have not examined it in detail, we suspect that F tc values from parametric XSPEC-tjpe 
fitting will also be inaccurate and biased except for very strong sources. 

We derive statistical and systematic errors on intrinsic source fluxes following the pro- 
cedure presented in §2.4.21 Figures [T5"b.fc compare the inferred statistical and systematic 
errors on total-band intrinsic source flux. Here, systematic errors based on X-ray model un- 
certainty generally exceed statistical errors; again note the expanded vertical scale in panel 
b. A maximal accuracy of ±60% (or ±0.2 dex in log flux) is achieved only for sources with 
> 50 total-band net counts. 

In contrast to total-band absorption-corrected fluxes, the systematic errors of hard- 
band intrinsic fluxes do not exceed statistical errors for < 50 hard-band count sources. Here, 
systematic and statistical errors become comparable for brighter sources. An accuracy of 
~ 60% for hard-band intrinsic fluxes can be achieved for sources with as few as ~ 7 — 10 hard- 
band net counts, unless very high absorptions (MedE t > 4 keV) are present. We thus find, 
as expected from qualitative considerations, that hard-band absorption-corrected intrinsic 



12 Considering that the log Nh values inferred nonparametrically using two-temperature simulated mod- 
els are systematically lower than that of the one-temperature XSPEC fitting ( §3.1|) . it is clear that one- 
temperature models can not fully recover the soft intrinsic PMS X-ray emission subject to absorption. 
This is further indicated by a comparison of 2-temperature m odel fits of the deep COUP dataset with 1- 
temperature model fits of an earlier 80 ks Chandra exposure (jFeigelson et al.l 12002). The COUP analysis 
gave ~ 2 times higher fluxes than the earlier analysis. This result emphasizes the dangers of absorption 
correction in fluxes that include the soft band. 
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fluxes are much more stable than intrinsic full-band fluxes. We thus encourage use of Fh c 
when absorption-corrected fluxes and luminosities are sought using nonparametric methods. 



4. Applications to Other Classes of Faint X-ray Sources 



Recall from §2.11 that a bright calibration sample of young stars in the Orion Nebula 
indicated a generic spectral model for this class of X-ray sources: a two-temperature thermal 
plasma with the cooler temperature fixed and the hotter temperature scaled to X-ray lumi- 
nosity (Figures [1] and [2]). This property of PMS stars allowed us to simulate X-ray models, 
and calibrate the derivation of astrophysical absorption and fluxes from the observed counts. 
Other classes of X-ray sources will exhibit different intrinsic spectra shapes with somewhat 
different dependencies on observable quantities. 

Spectral model families can be constructed for other classes of X-ray sources, extending 
nonparametric estimation methods beyond young star populations. Close binary star sys- 
tems with accreting black holes typically show two spectral components, a thermal disk and 
a nonthermal powerlaw. The relative strength of the compon ents scales with the source lumi- 
nosity and spectral hardness (jRemillard fc McClintockll2006l ). These relationships, based on 
brigh t Galactic X-ray binary systems, can be applied to X-ray b inary systems in nearby spiral 



(e.g. ^Griffiths et al 



2003 



Sarazin et al 



2000l : IPietsch et al.ll2005l : IStiele et al.ll2008f ) or elliptical flPrestwich et al. 



2003h galaxies. While the ultraluminous X-ray binaries in these galaxies 



may be bright enough for detailed parametric spectral modeling, the more typical source is 
well-adapted to our methods based on broad band counts and median energies. Spectral 
model families might also be con structed for faint X-ray source populations in the Galactic 
Center region (IMuno et al.ll2003l ). 



The X-ray spectra of faint active galactic nuclei (AGN) are often parametrized by a 
power law with index T subject to differing amounts of absorption, although more complex 
emission and absorption components may be present. At low redshift and luminosities, no 



statistical relationships are seen 



(jCappi et all 120061 ; IWinter et al. 



jet we en T, luminosity, and absorption in the X-ray band 



2008 ). But at higher redshifts, q uasar spectra appear to 



soften as luminosity increases (IShemmer et al.ll2006l ; ISaez et al.ll2008l ). As there is substantial 
scatter in T at both low and high redshifts, care must be taken to propagate uncertainties 
in spectral index through our procedures for estimating errors in the inferred broad-band 
fluxes and absorption. 
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5. CONCLUSIONS 

In this work, we show that important properties of X-ray sources — such as line-of- 
sight absorption, apparent broad-band fluxes, and intrinsic fluxes — can be estimated with 
reasonable accuracy using the easily measured photometric quantities X-ray count rate and 
median energy. The translation from observed to intrinsic quantities is achieved without 
employing often-uncertain and time consuming methods of non-linear parametric spectral 
modeling. While parametric modeling is best for strong sources where spectral shape details 
are well-populated with photons, we believe nonparametric estimates are preferred for faint 
sources. The concept of our procedures is similar to long-standing methods based on color- 
magnitude diagrams in optical and infrared astronomy. 

Specifically, we first demonstrate that the column density logiV# can be estimated 
directly from the total-band median energy MedE t with quantifiable uncertainties. The 
apparent total-band (0.5 — 8 keV) and hard-band (2 — 8 keV) fluxes can be estimated from 
the source count numbers NC t and NCh, and the absorption-corrected intrinsic fluxes F hc 
and F tc can then be derived. The methods are intended for sources with greater than ~ 5 — 7 
counts but fewer than 100 — 300 counts where parametric spectral fitting methods will be 
superior. 

We establish limitations to the method which are probably generally applicable for 
any reasonable spectral model family. Nonparametric procedures provide poor measures 
of absorption for sources with MedE t < 1.7 keV (or equivalently, log Njj < 22.0 cm -2 ). 
Absorption corrections for total-band fluxes, F tc , are often large and unreliable for sources 
with < 50 total-band counts. However, absorption corrections for hard band fluxes, Fh c , are 
reliable down to 7 - 10 hard-band countfi F hc estimation using nonparametric techniques 
has small biases and uncertainties only somewhat larger than optimal y/~N errors. 

Quantitative results on biases, statistical errors, and systematic errors for each estimator 
are given in the text and figures. These detailed results are linked to assumptions of the 
intrinsic spectral models of the source populations, which we treat here to be absorbed 



13 In the current implementation of the method the hard-band apparent flux (and thus hard-band net 
counts) as well as the full-band median energy are used as the primary input photometric quantities for 
derivation of intrinsic source properties (see Step 2 of our procedure described in 32.4.1[) . As discussed in 
the paper, several net counts in the hard-band are needed in order to derive meaningful intrinsic properties. 
This level of hard-band signal may be difficult to detect from an intrinsically soft X-ray source. Researchers 
interested in non-parametric characterization of intrinsically soft weak X-ray sources might consider modify- 
ing the current implementation of the method to use the total-band apparent flux (instead of the hard-band 
flux) as a primary input photometric quantity. 
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two-temperature thermal plasmas associated with PMS stars in star-forming regions. We 
calibrated the spectral models to the well-studied COUP sample, and successfully applied 
the methods to a new sample in the M17 star-forming region. 

The IDL implementation of our nonparametric method, XPHOT.pro, is provided. Com- 
puting time needed to derive intrinsic properties for ~ 2000 X-ray sources is < 1 minute. 
Our method and results can be directly used in statistical X-ray studies of young rich stel- 
lar populations; for example, for constructing column density maps and X-ray luminosity 
functions of young stellar clusters. Nonparametric results also can be used as initial and/or 
frozen parameters to use with the parametric XSPEC method for fine-tuning of spectral 
properties. 

The procedures described here are developed in two specific contexts, but are more 
broadly applicable. First, they use data products gener ated by the AC IS Extra ct data re- 



duction package designed for the Chandra ACIS detector (IBroos et al.ll2002l . l2009l ). However, 
the methods are not restricted to this software package or detector; the necessary quantities 
such as median energy and count numbers can be readily obtained using CIAO, MIDAS, 
FTOOLS, or other packages. The methods could be applied to X-ray CCD observations 
with the ASCA, XMM-Newton, Suzaku, or other space-borne X-ray telescopes if the results 
are re-calibrated to each mirror-detector combination. 

Second, the analysis is conducted for specific two-temperature thermal plasma spectral 
models associated with PMS stars, and our quantitative results are restricted to this appli- 
cation. However, similar analysis, likely resulting in similar results, can be made for other 
classes of sources, such as active galactic nuclei and X-ray binary star systems in nearby 
galaxies. Researchers extending these methods to such classes should explicitly define and 
simulate their own models to construct calibration tables to use with XPHOT.pro. 

We thank Michael Eracleous (Penn State) for helpful discussions, and the anonymous 
referee for helpful comments. This work is supported by the Chandra ACIS Team (G. 
Garmire, PI) through the SAO grant SV4-74018 and NASA Astrophysics Data Program 
grant NNX09AC74G (E. Feigelson, PI). 
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Table 1. Template X-ray Spectral Models for COUP stars 



log(ifcc) 


fcTi 


kT 2 




(crgss- 1 ) 


(kcV) 


(kcV) 




(1) 


(2) 


(3) 


(4) 


27.0 


0.8 


0.8 


0.2 


27.5 


0.8 


1.1 


0.4 


28.0 


0.8 


1.5 


0.6 


28.5 


0.8 


1.8 


1.0 


29.0 


0.8 


2.2 


1.3 


29.5 


0.8 


2.5 


1.9 


30.0 


0.8 


2.8 


2.5 


30.5 


0.8 


3.2 


2.8 


31.0 


0.8 


3.5 


3.2 


31.5 


0.8 


3.9 


3.2 


32.0 


0.8 


4.2 


3.2 



to 
o 



Note. — Column 1: Hard-band intrinsic lu- 
minosity. Column 2: Temperature of the cool 
plasma component. Column 3: Temperature of 
the hot plasma component. Column 4: Ratio of 
emission measures between hot and cool compo- 
nents. 



Table 2. Simulated Data to Figure [3] 



MedEt 
(kcV) 
(1) 


F phot /F 3 
(2) 


(ergs s" 1 ) 
(3) 


Band 

(4) 


0.971 


1.049 


27.0 


full 


0.986 


1.072 


27.0 


full 


1.000 


1.104 


27.0 


full 


1.029 


1.146 


27.0 


full 


1.059 


1.187 


27.0 


full 


1.102 


1.229 


27.0 


full 


1.146 


1.266 


27.0 


full 


1.190 


1.298 


27.0 


full 


1.234 


1.323 


27.0 


full 


1.292 


1.342 


27.0 


full 


Note. — 


Column 1: 


Apparent full-band median 



to 



energy. Column 2: Ratio of the apparent X-ray 
flux of simulated data from ACIS Extract photome- 
try to that inferred from flux integration of the in- 
put models evaluated on the energy bins set by the 
Chandra response matrix. Column 3: Hard-band 
intrinsic luminosity indicates corresponding simu- 
lated X-ray spectral model family from Table ^ 
Column 4: Energy band, "hard" or l 'full" . This ta- 
ble is available in its entirety in a machine-readable 
form in the online journal. A portion is shown here 
for guidance regarding its form and content. 



Table 3. Simulated Data to Figure H] 



MedE t 


log(JV H ) 




(keV) 


(cm- 2 ) 


(orgs s X ) 


(1) 


(2) 


(3) 


0.971 


20.000 


27.0 


0.986 


20.658 


27.0 


1.000 


20.977 


27.0 


1.029 


21.301 


27.0 


1.059 


21.477 


27.0 


1.102 


21.602 


27.0 


1.146 


21.699 


27.0 


1.190 


21.778 


27.0 


1.234 


21.845 


27.0 


1.292 


21.903 


27.0 



to 



Note. — Column 1 : Apparent 
full-band median energy. Column 2: 
X-ray column density. Column 3: 
Hard-band intrinsic luminosity indi- 
cates corresponding simulated X-ray 
spectral model family from Table ^ 
This table is available in its entirety in 
a machine-readable form in the online 
journal. A portion is shown here for 
guidance regarding its form and con- 
tent. 



Table 4. Simulated Data to Figure 



M edE t 
(kcV) 
(1) 


Log(iV-F) 
(2) 


log(ifec) 
(ergs s _1 ) 
(3) 


Band 

(4) 


0.971 


0.01311 


27.0 


full 


0.986 


0.06390 


27.0 


full 


1.000 


0.11212 


27.0 


full 


1.029 


0.23303 


27.0 


full 


1.059 


0.3301S 


27.0 


full 


1.102 


0.41713 


27.0 


full 


1.146 


0.49527 


27.0 


full 


1.190 


0.56576 


27.0 


full 


1.234 


0.62951 


27.0 


full 


1.292 


0.68732 


27.0 


full 



Note. — Column 1: Apparent full-band me- 
dian energy. Column 2: Intrinsic to apparent 
flux ratio of simulated data. These arc flux 
integrations of the input models (not rc-fittcd 
simulated spectra) evaluated on the energy bins 
set by the Chandra response matrix. Column 3: 
Hard-band intrinsic luminosity indicates corre- 
sponding simulated X-ray spectral model family 
from Table |T1 Column 4: Energy band, "hard" 
or "full". This table is available in its entirety 
in a machine- readable form in the online journal. 
A portion is shown here for guidance regarding 
its form and content. 



Table 5. Simulated Data to Figure [9] 



MedE 


AMedB 


l°g(£/ic) 


Band 


NC 


(kcV) 


(koV) 


(ergs s" 1 ) 




(cnts) 


(1) 


(2) 


(3) 


(4) 


(5) 


2.431 


0.2672 


27.0 


hard 


5 


2.445 


0.2716 


27.0 


hard 


5 


2.460 


0.2784 


27.0 


hard 


5 


2.475 


0.2868 


27.0 


hard 


5 


2.489 


0.2923 


27.0 


hard 


5 


2.504 


0.2981 


27.0 


hard 


5 


2.518 


0.3052 


27.0 


hard 


5 


2.533 


0.3131 


27.0 


hard 


5 


2.548 


0.3180 


27.0 


hard 


5 


2.562 


0.3212 


27.0 


hard 


5 



Note. — Column 1: Apparent median energy in the 
corresponding energy band. Column 2: Statistical er- 
ror on median energy. Column 3: Hard-band intrin- 
sic luminosity indicates corresponding simulated X-ray 
spectral model family from Table ^ Column 4: En- 
ergy band, "hard" or "full". Column 5: Net counts in 
the corresponding energy band. This table is available 
in its entirety in a machine-readable form in the online 
journal. A portion is shown here for guidance regarding 
its form and content. 
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Fig. 1. — Spectral input models (squares) for our simulations based on the COUP data 
(circles) of ~ 5 00 Orion Nebula Cl uster (ONC) PMS stars with available 2-temperature 
model fits from iGetman et al.l (120051 ). The hot temperature component (panel a) and ratio 
of emission measures between hot and cool components (panel b) are plotted against intrinsic 
hard-band (2.0 — 8.0 keV) X-ray luminosity. The dashed curves encompass ~ 75% of the 
ONC stars and are constructed such that at each log(L^ c ) bin their y-axis value is equal to 
the y-axis value of the "second nearest neighbor" simulated (square) model; that is, the top 
(bottom) dashed curves are versions of the solid curves shifted two log(L/ lc ) bins to the left 
(right). 
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Fig. 2. — Same as Figure [H but restricted to the CO UP PMS star sub-sample with an 



available accretion indicator (e.g. iPreibisch et al.ll2005l ). COUP stars classified as active 



accretors having the Ca II 8542 A line in emission with equivalent width EW(Ca II) < — 1A 
are shown with • symbols, and stars classified as weakly accreting or non-accreting with 
absorption equivalent width EW(Ca, II) > 1A are shown with x symbols. 
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1 2 3 4 5 6 

MedE (0.5-8keV band) [keV] 

Fig. 3. — Ratio of the apparent X-ray flux of simulated data from ACIS Extract photometry 
to that inferred from flux integration of the input spectral models plotted against X-ray 
median energy in the (a) total band (0.5 — 8.0 keV) and (b) hard band (2.0 — 8.0 keV). 
Results are presented for all 11 simulated spectral model families with the model families 
for \ogL hc = 27 ergs s" 1 ()ogL hc = 32 ergs s _1 ) stars as the top (bottom) curves. 
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1 2 3 4 5 6 

MedE (0.5-8keV) [keV] 



Fig. 4. — Simulated calibration curves of X-ray column density plotted against X-ray median 
energy, shown for the 11 simulated spectral model families from Figure [TJ The top (blue 
dashed) and bottom (blue solid) curves correspond to the model families for log(L/ lc ) = 
27 ergs s _1 and log(Lhc) = 32 ergs s _1 stars, respectively. The insert exemplifies estimation 
of the column density, \ogN HtPhot , and its statistical (red arrows) and systematic (green 
arrows) errors. In this example, the source has 20 counts with median energy 2.5 keV and 
median energy uncertainty 0.35 keV, using the log(L/ lc ) ~ 29 ergs s _1 spectral model family. 
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1 2 3 4 5 6 

MedE (0.5-8keV) [keV] 



Fig. 5. — Intrinsic to apparent flux ratio of simulated data are plotted against X-ray median 
energy. Results are shown for the 11 simulated spectral model families (solid lines) with 
\og(L hc ) = 27 ergs s" 1 and \og(L hc ) = 32 ergs s" 1 as the leftmost and rightmost curves, 
respectively. Model curves are truncated to indicate the locus of ~ 1400 ONC PMS stars 
from the sensitive COUP project. Correction factors are shown for the (a) total band 
(0.5 - 8.0 keV) and (b) hard band (2.0 - 8.0 keV). 
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Step 1 . Read Acis Extract 
photometry outputs. 



AE 




MedEt 
MedEh 
NCt 
NCh 




Step 2. Obtain "best model family" (BMF). 



Pass through Equation 4. 
Get apparent hard-band 
luminosity (so far slightly 
biased). 




(LhjDiase^) 



Iterate through Fhc/Fh 
versus MedEt simulated 
data (Table 4) to get BMF. 



BMF - Lh_biased 



Step 3. Debias photometric fluxes; get 
apparent fluxes in units of ergs/cm 2 /s; 
get absorbing column density. 




Pass through 
Equations 3-4 
and Flux Bias 
versus MedEt 
simulated 
data (Table 2). 




MedEf 
vBMF 



Pass through 
N u versus MedEt 

rl 

simulated 
data (Table 3). 




Step 4. Obtain intrinsic fluxes 
in units of ergs/cm 2 /s. 

fh~ 
Ft 

MedEt, 
sBMF 



Pass through intrinsic to apparent 
flux ratio versus MedEt simulated 
data (Table 4). 




Fig. 6. — The procedure of our nonparametric method given in the form of a flowchart. The 
procedure for obtaining intrinsic fluxes and absorbing column density. 
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Step 5. Obtain statistical errors on 






/NCh \ 
/ NCt > 


median energies and apparent fluxes. 




MedEt 




Pass through AMedE versus MedE 






\ MedEh J 


— > 


simulated data (Table 5) to get 








errors on median energy. 




/NCh \ 
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MedEh 
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\ aNCt / 
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— > 


Error propagation in Equations 
1 -4 to get statistical errors on 
apparent fluxes. 


> / aFt\ 
\aFh/ 












Step 6. Obtain statistical and systematic errors 
on absorbing column density and intrinsic fluxes. 



Pass through N H versus MedEt simulated 
data (Table 3) to obtain errors on 
absorbing column density. 




Pass through intrinsic to apparent 
flux ratio versus MedEt simulated 
data (Table 4) to get errors on 
intrinsic fluxes. 



*jFhc,staf\ 
aFhc,syst 
aFtcstat 
fFtcsyst, 



Fig. 7. — The procedure for obtaining errors on intrinsic fluxes and absorbing column density. 
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Fig. 8. — X-ray "color- magnitude diagram" derived from our method with apparent hard- 
band flux versus apparent total-band median energy. Lines represent X-ray model isochrones 
derived from our nonparametric method: dashed lines indicate isochrones of intrinsic hard- 
band X-ray luminosity, solid lines indicate isochrones of absorbing column density. 
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Fig. 9. — Uncertainties of the estimated X-ray median energy plotted against median energy 
for different net count strata from simulated data in the (a) total band (0.5 — 8.0 keV) and 
(b) hard band (2.0 — 8.0 keV). For each count stratum, shown with a different color, results 
are shown for the 11 simulated spectral model families from Figure [T]with log(L/ lc ) = 27 ergs 
s _1 and log(L? iC ) = 32 ergs s _1 as the bottom and top curves, respectively. Moving along 
an individual model family curve from left to right corresponds to increasing the column 
density parameter of that spectral model family (see Figure 0J. 
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Fig. 10. — (a) X-ray column density inferred from XSPEC spectral 1-temperature model 
fits of > 1600 M17 sources plotted against X-ray median energy. Two solid curves are 
calibration lines corresponding to the simulated model families of log(L/ lc ) = 29 ergs s _1 
(top) and \og(L hc ) = 31 ergs s" 1 (bottom) from Figure HI (b) Ratio of X-ray column density 
obtained from nonparametric estimation to that inferred from XSPEC spectral fits of the 
Ml 7 data plotted against X-ray median energy. Net count strata are color-coded. 



1.1 

1.0 
0.9 
0.8 



x 7 

CO 



1(a) 



i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i 







a 0-6 h 
i 0.5 

b 9 0.4 
0.3 
0.2 
0.1 
0.0 
1.1 
1.0 
0.9 
0.8 

1 0J 
^1 0.6 h 

CO 

| 0.5 
I 0.4 

0.3 

0.2 

0.1 



o m 



i i i i i i i i . i 




®fc^"S7T" " r • " • °i 

Till" I I I . 



1(b) 



0.0 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



NC=5-7 

7-10 

10-20 

20-50 

50-100 

100-300 

>300 



I I t t I I . I J. i U I I ii I 



I I I 







2 3 4 

MedE (0.5-8.0 keV) [keV] 



Fig. 11. — Statistical (a) and systematic (b) errors on the column density log Nh inferred 
nonparametrically plotted against X-ray median energy. Net count strata are color-coded. 
Results are presented for ~ 2000 M17 sources. 




Fig. 12. — Ratio of apparent (not corrected for absorption) de-biased photometric X-ray 
flux to that inferred from XSPEC spectral fitting of Ml 7 data plotted against X-ray median 
energy. Results are presented for ~ 1600 M17 sources for the total band (panel a), and 
~ 1300 M17 sources for the hard band (panel b). Net count strata are color-coded. Fluxes 
are in units of ergs s _1 cm -2 . 
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Fig. 13. — Statistical errors on (a) the total band and (b) the hard band apparent X-ray 
fluxes inferred nonparametrically plotted against X-ray median energy for ~ 2000 and ~ 1800 
Ml 7 sources, respectively. Net count strata are color-coded. 
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Fig. 14. — Comparison of intrinsic fluxes (absorption-corrected) obtained from our nonpara- 
metric methods with those inferred from XSPEC spectral fits (panels a, c), and their ratios 
against X-ray median energy (panels b, d). Results are presented for ~ 1600 M17 sources 
for the total band (upper panels), and ~ 1300 M17 sources for the hard band (lower panels). 
Net count strata are color-coded. 
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Fig. 15. — Statistical (left panels) and systematic (right panels) errors of nonparametrically 
derived intrinsic fluxes corrected for absorption. Results are presented for ~ 2000 M17 
sources for the total band flux (upper panels) and for ~ 1800 M17 sources for the hard band 
(lower panels). Net count strata are color-coded. To avoid symbol clutter only 5 — 7 and 
50 — 100 count strata are presented on panels (b) and (d). 



